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ABSTRACT 

We present deep XMM-Newton EPIC spectroscopic and timing X-ray observations of the nearby 
solitary radio millisecond pulsar, PSR J0030-I-0451. Its emission spectrum in the 0.1-10 keV range 
is found to be remarkably similar to that of the nearest and best studied millisecond pulsar, PSR 
J0437~4715, being well described by a predominantly thermal two-temperature model plus a faint 
hard tail evident above ^2 keV. The pulsed emission in the 0.3-2 keV band is characterized by two 
broad pulses with pulsed fraction ^--^60-70%, consistent with a mostly thermal origin of the X-rays only 
if the surface polar cap radiation is from a light-element atmosphere. Modeling of the thermal pulses 
permits us to place constraints on the neutron star radius of i? > 10.7 (95% confidence) and R > 10.4 
km (at 99.9% confidence) for M = 1.4 Mq. 

Subject headings: pulsars: general — pulsars: individual (PSR J0030-)-0451) — stars: neutron — 
X-rays: stars — relativity 



1. INTRODUCTION 

PSR J0030-I-0451 is one of the nearest known rotation- 
powered "recycled" millisecond pulsars (MSPs) in the 
field of the Galaxy {D = 300 ±90 pc; Lommen et al. 
2006), with a spin period P = 4.87 ms and intrinsic 
spindown rate P = dP/dt = 1.0 x 10~^° s s~-^, implying 
a surface dipole magnetic field strength B « 2.7 x 10® 
Gauss, a characteristic age r « 7.8 Gyr, and spin-down 
luminosity i? w 3 x 10^^ ergs s^^ This sohtary MSP 
was discover ed at radio frequenci es in the Arecibo drift 
scan survey (jLommen et al.|[2000l ). Recently, a firm de- 
tection of PSR J0030-I-0451 at 7-ray ener gies by the 
Fermi Large Area Telescope was reported bv lAbdo et al.l 
(|2009f) . In X-rays, PS R J0030-I-0451 has b een detected 
with R OSAT PSPC (iBecker et al.l [200I and XMM- 
Newton ijBecker fc AschenbachI 120*021 ) . Despite the lim- 
ited photon statistics, these observations clearly showed 
two distinct broad pulses. In addition, the XMM-Newton 
observation revealed a relatively soft 0.3-2 keV spe ctrum, 
quali tatively similar to that of PSR J0437-4715 (jZavlinl 

X-ray observations have detected a number of rotation- 
powered MSPs that exhibi t soft, presumably thermal, 
emission (jZavhni i2006; B ogdanov et al.l l2006at IZavlinI 
l2007f ). This radiation likely originates from the pul- 
sar magnetic polar caps that are heated by ener- 
getic particles from the pu lsar magnetosphere (see, e.g., 
iHarding fc Muslimovll2063 ). As this heat is confined to 
a small portion of the NS, study of the X-ray proper- 
ties of MSPs could offer insight into key NS proper- 
ties that are inaccessible by other observational means 
(e.g., radio pulse timing) such as the radiative prop- 
erties of the NS surface, magnetic field geometry, and 
NS compactness {R/Rs, where Rs — 2GM/c^ and R 
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an d M are the stellar mass and radius). As shown 
bvlPavlov fc Zavlin (1997^■ IZavhn fc Pavlov! (|199^ . and 
iBogdanov et al. (2007 ), a model of polar cap thermal 
emission from an optically-thick hydrogen (H) atmo- 
sphere provides a good description of the X-ray pulse 
profiles of PSR J0437~4715, the nearest known MSP. 
On the other hand, a blackbody model is inconsistent 
with the pulsed X-ray emission. Furthermore, there 
is compelling evidence for a magnetic dipole axis off- 
set from the NS center ()Bogdanov et al.ll2007l) . Finally, 
the compactness of PSR J0437-4715 is constrained to 
be R/Rs > 1.6 (99.9% confidence), whi ch for the cur- 
rent b est mass measurement fl.7 6 Mr^. see|Verbe ist et all 
I2008f ) implies R > 8.3 km. Bog danov et all (|2008l have 
shown, using short archival XMM-Newton observations 
of PSRs J0030+0451 and J2124-3358, that even with rel- 
atively crude photon statistics, it is possible to set limits 
on M/R. In particular, assuming 1.4 Mq, the stellar ra- 
dius is constrained to be i? > 9.4 km and R > 7.8 km 
(68% confidence) for each pulsar, respectively. Thus, re- 
alistic modeling of X-ray data from MSPs appears to be 
a promising approach towards revealing the structure of 
NSs and warrants further X-ray observations of nearby 
MSPs. As one of the nearest known recycled pulsars, 
PSR J0030-I-0451 is well suited for such an investigation. 

In this paper, we present an analysis of deep XMM- 
Newton spectroscopic and timing observations of PSR 
J0030-I-0451. These observations provide new insight 
into the X-ray properties of MSPs and the structure of 
neutron stars. The work is presented as follows. In §2 
we summarize the observations and data reduction pro- 
cedure. In §3 we briefly examine the X-ray images. In §4 
we constrain the spectral properties of the pulsar, while 
in §5 we conduct a timing analysis. In §6 we attempt 
to model the pulsed emission from the MSP. We offer 
conclusions in §7. 

2. OBSERVATION AND DATA REDUCTION 

PSR J0030-f0451 was targeted by XMM-Newton be- 
tween 12 and 14 December 2007 (observation ID 050229). 
The observation was carried out in a single uninterrupted 
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Fig. 1.— XMM-Newton MOSl and MOS2 mosaic image of the 
field near PSR J0030+0451 in the 0.3-2 keV band. The rectangle 
shows the relative position and orientation of the active EPIC pn 
chip used in this observation, with the timing redout direction in- 
dicated by the arrow. The dashed lines correspond to the EPIC pn 
detector columns used to extract the source counts from the pulsar 
while the circles mark the positions of sources that fall within the 
EPIC pn field of view. 

130-kilosccond exposure, corresponding to the entire us- 
able observing time of XMM-Newton revolution 1467. 
The European Photon Imaging Camera (EPIC) MOSl 
and M0S2 instruments were configured for full imaging 
mode. The EPIC pn was configured for fast timing mode, 
in which only CCD 4 is active, allowing 30 /iS time reso- 
lution at the expense of one imaging dimension. For all 
three instruments, the thin optical filter was used. Due 
to the faint nature of J0030+0451, the dispersed Reflec- 
tion Grating Spectrometer (RGS) data provide no useful 
spectral or timing information. Therefore, the grating 
data were not used in our analysis. 

The data reduction, imaging, and timing analyses were 
performed using the SAS'' 8.0.1 and FTOOLS^ 6.6.1 soft- 
ware packages, while the spectral analysis was conducted 
in XSPEC^ 12.5.0. The raw MOS and pn datasets were 
first reprocessed using the SAS emchain and epchain 
pipelines, respectively, and subsequently screened for in- 
stances of high background proton flares. A total of ^^36 
ks of high-background data were discarded, resulting in 
92.5, 92.1, and 92.1 ks of usable exposure time for MOSl, 
M0S2, and pn, respectively. Finally, the appropriate 
pattern, pulse invariant, and flag filters were applied to 
produce data suitable for spectral and timing analyses. 

The XMM-Newton SAS is developed and maintained by the 
Science Operations Centre at the European Space Astronomy Cen- 
tre and the Survey Science Centre at the University of Leicester. 

^ Available at http://heasarc.gsfc.nasa.gov/ftools/ 

^ Available at http://heasarc.nasa.gov/docs/xanadu7xspec7] 




20 40 
Detector column 



Fig. 2. — The distribution of counts from XMM-Newton EPIC 
pn along the imaging dimension of the fast timing mode in the 0.5— 
2 keV band. The pair of dashed lines mark the detector columns 
that the pulsar counts were extracted from. The labels show the 
positions of the sources from Figure 1. 

For the spectroscopic analysis, source counts in the 
MOS 1 and 2 data were extracted from circles of radius 
60" centered on the radio MSP position, which enclose 
~90% of the total energy at 1.5 keV. The background 
was taken from three source- free regions surrounding the 
pulsar. For the spectroscopic and timing analyses of the 
EPIC pn dataset, the source counts were obtained from 
RAWX detector columns 31-41 (inclusive), equivalent to 
a radius of 22.5" in the RAWX direction, which encloses 
^74% of the total energy for 0.3-2 keV. This relatively 
narrow region was chosen in order to avoid contamina- 
tion from sources A and B as well as to minimize the 
large background level, which dominates beyond ^--^20" 
from the source position, arising due to the 1-dimensional 
imaging mode used for the pn detector. To ensure a 
reliable estimate, the background was taken from three 
source- free regions. 

3. IMAGING ANALYSIS 

Figure 1 shows the 0.3-2 keV coadded image from 
EPIC MOSl and M0S2 of the field in the vicinity of 
PSR J0030-f0451, as well as the relative orientation and 
readout direction of the EPIC pn detector. It is clear 
that the pulsar is the brightest source in the image. Due 
to the favorable telescope roll angle during the observa- 
tion, the EPIC pn detector columns containing the bulk 
of counts from the pulsar are effectively free of contam- 
ination from other sources (Fig. 2). Sources A and B 
combined contribute with only a few counts (<0.1% of 
the total) in the source extraction colums of the pulsar. 

In the MOSl and 2 images, the source count distri- 
bution around the radio position of the pulsar is fully 
consistent with that of a point source. Thus, we find no 
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indication of any diffuse extended emission that could 
arise due to a bow shock, as seen in PSR J2124-3358 
(|Hui fc Beckeij[2"006[ ). This is not surprising given that 
for the combination of spin-down luminosity (i? « 3 x 
10^'^ ergs s~^) and remark ably low transverse s pace ve- 
locity (8 — 17 km s~^, see lLommen et aT1l2006( ) for this 
pulsar, no detecta ble diffuse X-ray emission i s expected 
(see, e.g., Fig. 1 in lChatteriee fc Cordesll2002D . 

4. SPECTROSCOPIC ANALYSIS 

To facilitate the spectral fitting of the MOSl and 
M0S2 data, the extracted counts in the 0.1-10 keV band 
were grouped so as to ensure at least 30 counts per bin. 
For the pn data the counts in the 0.3-3 keV were grouped 
with at least 150 counts per bin. Due to the overwhealm- 
ing background level, no useful spectral information is 
present above 3 keV in the pn data. In addition, pho- 
tons in the 0.42-0.5 keV range were not used in order 
to eliminate an instrumental noise artifact specific to the 
fast timing mode of the pn detector. 

As foun d in previo u s stud ie s of PSR J003 -I-0451 
jBcckcr fc AschenbachI l2002t iBogdanov eFaD 120081 : 
[Zavlin ,200^ , the continuum X-ray emission from this 
pulsar cannot be adequately described by a single com- 
ponent emission model. Consequently, in the spectral 
fits we apply the following plausible multi-component 
models: i) two thermal, ii) a single thermal plus 
non-thermal, iii) two Comptonized thermal iv) two 
thermal plus one non-thermal, and v) three thermal 
components. For the thermal components we consider 
both a blackbody model and the unma gnetizcd (B — 0] 
H atm osphere model first presented in iMcClintock et al l 
()2004f l and included as a private model in XSPEC. The 
assumption of _B = is appropriate since for typical 
MSP magnetic fields (~10^~^ G) the magnetic field does 
not affect the radiativ e opacities and equation of state 
of the atmosphere (see lRomanilll987L FZavlin et al.lll996l : 
IMcClintock et al.ll2004 for details). 

It has been found that blackbody fits to MSP spec- 
tra result in effective emission areas that are an order 
of magnitude smaller than the expected polar cap area, 
Rpc = {2iTR/cPy^^R, whil e for a H atmosphere they 
are found to be comparable (iBecker fc Aschenbachll2002l: 
iZavlin et al.ll2002l: IBogdanov et al.l l2006a: Zavlin 2006). 
The discrepancy could, in principle, arise due to non- 
uniform heating of the polar caps. More plausably, it 
suggests that the surface of the polar caps is covered by 
a light-element atmosphere. As the spectrum of a H at- 
mosphere model is harder than that of a blackbody for 
the same effective temperature, a blackbody model fit- 
ted to a H atmosphere continuum would yield a higher 
inferred temperature and, as a result, a smaller effective 
area. Nonetheless, for the sake of completeness and eas- 
ier comparison with previous X-ray studies of MSPs, we 
have considered blackbody models in our spectral fits as 
well. 

IZavlin et al.l ([2002), 
IBogdanov et al.l p007. ) have noted that due to the 
energy-dependent limb darkening effect for the atmo- 
sphere model, even for phase-integrated spectra it is 
necessary to take into account the rotation of the star as 
well as the relative orientation of the spin axis, hot spots, 
and the line of sight to the observer in order to ensure 
a reliable measurement of the temperature and emission 
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Fig. 3.— The XMM-Newton EPIC X-ray spectra of PSR 
J0030-I-0451 fitted with a two-temperature H-atmosphere plus 
powerlaw model. The middle panel shows the fit residuals while 
the bottom panel shows the best fit model. See Table 1 for best fit 
parameters. 

radiu s (see, in particular, Fig. 2 of IBogdanov et al.l 
[200l. For this reason, we employ the spectral model 
from IBogdanov et al.l (|2007| ). This model allows us to 
specify the angles between the spin and magnetic axes 
[a) and the spin axis and the line of sight (C) as well as 
the mass and radius of the pulsar in the spectral fits. 
Thus, the normalization of the thermal component(s) 
yields the true (deprojected and redshift-corrected) 
effective area of the emission region(s). In addition, the 
quoted best-fit temperatures represent redshift-corrected 
values, i.e. as they would be measured at the NS surface. 
For the calculation of all effective radii, we consider 
the parallax distance D — 300 pc. In the case of the 
H atmosphere model, we assume a neutron star with 
M = 1.4 M© and i? = 12 km. Based on the X-ray 
pulse profile morphology (discussed in §5) we also fix 
a — 70° and C, — 80° since this combination of angles 
can reproduce the X-ray pulse profile shape for M = 1.4 
M© and i? = 12 km. All uncertainties quoted below are 
given at a la level for one interesting parameter. The 
various spectral models and their best fit parameters 
are summarized in Table 1. 

4.1. Two- Temperature Thermal Spectrum 

First, we consider the two-temperature thermal 
model that was found to provide a good fit to 
the archival XMM-Newton d a ta of PSR J003 0-I-0451 
([Becker fc AschenbachI 120021: IBogdanov et~an [2008.) . 
This model describes the new spectrum well up to ~3 
keV but cannot fully account for the excess flux at higher 
energies. Thus, it is necessary to invoke alternative emis- 
sion models. 

4.2. Thermal Plus Non-thermal Spectrum 

Next, we apply a composite thermal plus non-thermal 
(powerlaw) model. As evident from Table 1, for either 
a blackbody or H atmosphere thermal component the 
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TABLE 1 

Best fit spectral models and unabsorbed fluxes for PSR J0030+0451. 



Model^' 


(1020 cm-2) 


Tcft.l 
(lO'' K) 


fleff.l*^ 
(km) 


^"011,2 

(106 K) 


(km) 






Fx" 
(0.1-10 keV) 




BB(x2) 
Hatm(x2) 


0.90+0 3 

2 0+°-'^ 
'^■"-0.3 


-, 71 +0.16 


'J-'JO-0.03 
'-'■^^-0.03 




9 9+0.7 
^•^-0.6 






3.1 ±0.1 
3.7 ±0.1 


1.21/299 
1.16/299 


r 


BB+PL 

Hatm+PL 

BB(x2)+PL 

Hatm(x2)+PL 


4+0.4 
'-'■^-0.4 

n O + 0.5 

3.iS_0 7 

1 5+°-" 
1 Q+0.7 

j^-y-0.6 

1 S+0-2 


9 1 1 +0.U5 
^■^^-0.04 

1 06+°°^ 

2 04+0.86 

2 P-,+0.11 
^•^^-0.11 

1 CI +0.08 
'—^'--0.07 


o.orl^i 

0.08lH^ 
n 07+0.02 
'-'•■^'-0.01 
3+°-2 


, yc+o.aa 

^- '"^-0.24 
, -,0+0.08 
J^-^°-0.08 

U.D9_o 08 

n 70+0-08 

"J-0.07 


^•J^J^-0.13 

•^A+O-ao 

^•■^^-0.17 

1 n+0.19 
1-^-0.16 


■^•'-'^-0.08 
06+0-09 

ns+0-09 

1.4 

2 n+0.2 

1.4 




8.6 ±0.1(7.1 ±0.1) 
6.9 ± 0.1(4.8 ± 0.1) 
8.2 ± 0.2(6.6 ± 0.2) 

3.7 ±0.2(0.5 ±0.1) 
3.9 ± 0.4(0.6 ± 0.1) 
3.7 ±0.3(0.4 ±0.1) 


1.12/299 
1.11/299 
1.10/297 
1.13/298 
1.12/297 
1.11/298 














kTe 

(kcV) 


T 






Compbb(x2) 


-^•«-0.4 


9 90+U.iU 
^■^"-0.10 


'J--'- -"--0.05 




0.38±li;?^ 


100 


0.12 ±0.01 


3.8 ±0.3 


1.10/298 














Tcti,3 

(IQO K) 


^?cff.3 
(km) 






BB(x3) 
Hatm(x3) 


^■^-0.5 
9 9+0.6 


og+0:84 


u.uul_o 003 

03+° ''^ 
*J-^'^-0.03 


9 97+0.13 
^•^'-0.14 
1 05+0-07 

^■^'^-o.og 


0.95«;* 


^•'-'^-0.14 

40+0-06 


'^•™-0.29 

4 q+2-8 


3.7 ±0.3 
4.0 ±0.4 


1.10/297 
1.11/297 



^ PL is a powerlaw, BB a blackbody, and Hatm a H atmosphere model. All uncertainties quoted are la}" RcB calculated assuming 
a distance of 300 pc. For the H atmosphere model, the numbers quoted represent deprojected and redshift-corrected effective 
radii of one polar cap of a neutron star with AI = 1.4 Mq, R = 12 km, a = 70°, and ^ = 80° (see text for details).'^ Unabsorbed 
X-ray flux (0.1-10 keV) in units of 10"^'' ergs cm"^ s"^. The value in parentheses represents the flux contribution of the 
powerlaw component. 



model yields statistically acceptable fits. However, the 
derived values of iVn are significantly larger than those 
inferred from independent measurements. In particular, 
the total H column density through the Galaxy along 
the lin e of sight to the pulsar is only Nn « 3 x 10^*^ 
cm~^ (jPicke fc Lockmanlll99Cl[) . On these grounds, we 
deem the validity of the single temparature thermal plus 
non-thermal model doubtful. 

4.3. Comptonized Thermal Spectrum 

iBogdanov et al.l ()2006 b') have proposed that the hard 
tail seen in the X-ray spectrum of PSR J0437-4715 may 
be produced by inverse Compton scattering (ICS) of the 
soft thermal X-rays by relativistic of low optical depth 
(r <C 1). Such particles are likely present in the pul- 
sar magnetosphere above the polar caps. In the spectral 
fits of PSR J0030-I-0451, we employ the compbb Comp- 
tonized blackbody model (jNishimura et al.lll986f) . which 
includes two additional parameters: a scattering particle 
temperature kT^ and optical depth r. The model as- 
sumes a plane-parallel, semi-infinite scattering medium 
with the source of thermal photons at the bottom. The 
formalism employed is valid between kT < kT^ < 150 
keV. For the two thermal components we link the r pa- 
rameter as both would likely be scattered by the same 
population of particles. For a choice of kTe = 100 
kcV, the best fit parameters are Nn — 1.8to'4 x 10^*^ 

cm-2, Teff,! = {2.28tVw) x 10' K, R.s.i = O-lltaos 
km, reff,2 = (l-OStVol) X 10' R2 = OMt^oil km, 
r = 0.12 ± 0.01 with xl = 1-10 for 298 degrees of free- 
dom. The implied unabsorbed flux is 3.8 x 10"^'^ ergs 
cm~^ s~^ (0.1-10 keV). As is generally the case for ICS, 
kTe and r are strongly correlated so acceptible fits were 
obtained for a wide range of these parameters. Since 
we expect kT ^ kT^, the choice of fcTg as well as the 
actual energy distribution of the scattering particles do 



not significantly affect the X-ray spectral fits. Thus, the 
ICS model appears to provide a good description of the 
spectrum of PSR J0030-I-0451 as weU. 

4.4. Two- Temperature Thermal Plus Powerlaw 
Spectrum 

IZavlin et al] (|2002l ) andlZavlid (|2006l ) have found the 
continuum X-ray emission from PSR J0437-4715 to be 
well described by a two-temperature thermal and a single 
powerlaw components. As expected based on the similar 
spin properties of the two pulsars, the same model pro- 
vides a good description of the X-ray spectrum of PSR 
J0030+0451. Formally, an F-test indicates a 99.9997% 
probability that the addition of the powerlaw component 
to the two-temperature thermal model is required. As 
with the single temperature composite model, the de- 
rived best fit H column density for the two blackbody 
plus powerlaw model with F unconstrained are in excess 
of the upper limit A'h = 3 x 10^° cm~^ and yield a rather 
steep spectral photon index (F ~ 3). On the other hand, 
the two-temperature H atmosphere plus powerlaw model 
(Fig. 3) results in a more plausible value of Nn and a 
value of F very similar to that found for PSR J0437- 
4715 (F « 2). We note that due to the poor photon 
statistics above ~4 keV, statistically acceptable fits are 
obtained for a wide range of photon indices (1.0-3.5 at 
3cr confidence) although for F > 2.5 the implied Nh also 
exceeds 3 x 10^° cm~^ even for the two-temperature H 
atmosphere plus powerlaw model. 

4.5. Three- Temperature Thermal Spectrum 

The high-energy portion (>2 kcV) of the PSR 
J0030-f0451 X-ray spectrum could be due to an ad- 
ditional (third) thermal component. Indeed, for both 
three-temperature blackbody and H atmosphere models 
we obtain statistically acceptable fits. The best fit pa- 
rameters for this model suggest the presence of a very 
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Fig. 4. — Optical-to-7-ray spectrum of PSR J0030+0451 showing 
the XMM-Newto n X-ray data pres ented in this paper, VLT opti- 
cal upper hmits (Koptsevic h et al.l 120031) . and Fermi 7-ray mea- 
surements (Abdo et al. 200lD {green circles). The sohd blue line 
shows the total absorbed two-temperature H atmosphere plus pow- 
erlaw model spectrum for T = 2.0, while the red lines shows the 
best fit powerlaw model (solid) and uncertainties (dotted) for the 
Fermi LAT data. The dashed lines show the extrapolation of the 
best fit powerlaw spectrum with F = 2.0 and F = 1.4 of the two- 
temperature H-atmosphere plus powerlaw model to the optical and 
7-ray ranges. 

hot and very small (Tcff = 7 MK and R^s = 4 m for 
a blackbody or Teg = 3.4 MK and i?eff = 30 m for a 
H atmosphere) emission region on the neutron star sur- 
face. In principle, such a region could arise due to highly 
non-uniform magnctosphcric heating across the polar cap 
surface. 

4.6. Optical-to-gamma-ray Spectrum 

'Koptsev ich et al] (|2003D have conducted deep VLT 
F0RS2 BVR imaging of the field surrounding the ra- 
dio timing position of PSR J0030-I-0451. No plausible 
optical counterpart was detected down to B > 27.3, V 
> 27.0 and R > 27.0 in the i mmediate vicinity o f the 
radio pulsar position. In "/-ravs- lAbdo et al] (|2009( ) have 
reported on the Fermi LAT detection of pulsations from 
PSR J0030-f 0451 above 100 MeV. The 7-ray spectrum 
(100 MeV to 10 GeV) is well described by a powerlaw 
with a high-energy exponential cutoff with F = 1.4 ± 0.2 
and a cutoff energy of Ec = 1.7 ± 0.4 GeV. 

It is interesting to examine how the X-ray spectrum 
extrapolates to optical and 7-ray energies for the dif- 
ferent spectral models discussed above. Figure 4 shows 
the optical-to-7-ray spectrum of PSR J0030+0451 based 
on the available XMM-Newton, VLT, and Fermi obser- 
vations. The thermal components seen in the soft X- 
ray band contribute neglegibly to the optical and 7-ray 
fluxes. On the other hand, for the best fit two-component 
H atmosphere thermal plus power-law model, the spec- 
tral photon index of F = 2 grossly overestimates the 
optical fiux (by a factor of 30 — 40), violating the VLT 
upper limits, while it greatly underestimates the 7-ray 
flux (by a factor of ~100). This implies that such a pow- 
erlaw must break to a flatter powerlaw both in the UV 
range and between 10 keV and 100 MeV. 

Surprisingly, an extrapola tion of the best fit 7-ray spec- 
trum with F = 1.4 (from lAbdo et all l2009f ) is a very 
close match to the best fit two-compenent thermal plus 
non-thermal X-ray spectrum with F = 1.4 (see Table 



1). Within the uncertainties of the powerlaw normal- 
izations as well as the ~20% effective are a uncertainty 
of the Fermi LAT (see I Abdo et al.ll2009[ ). the two are 
consistent with being identical. This suggests that if the 
same emission process is responsible for the hard X-ray 
tail and the 7-ray emission from the pulsar, then the pow- 
erlaw index has to be very close to F = 1.4. However, 
given the rather large uncertainties in the two spectra, 
substantial improvement in the photon statistics in both 
energy ranges and detections at intermediate energies (10 
keV to 100 MeV) are required to offer more conclusive 
statements regarding the relation between the 7-ray and 
any non-thermal X-ray emission. 

Note that for the Comptonizcd thermal interpretation 
of the X-ray spectrum, one does not expect any contribu- 
tion of the model components in the UV/optical nor in 
the Fermi band as the high-energy cutoff of the Comp- 
tonizcd tail is likely well below the low energy threshold 
of Fermi LAT (~100 MeV). Thus, for this model there 
is no restriction on the slope of the hard tail. 

5. TIMING ANALYSIS 

To investigate the rotation-induced modulations of the 
X-ray flux from PSR J0300+0451, the photon arrival 
times extracted from the EPIC pn source region shown in 
Figure 2 were flrst translated to the solar system barycen- 
ter using the SAS barycen tool. The arrival times were 
subsequently folded at the pulsar spin period using the 
TEMPD^ pulsar timing p ackage and the rad io ephemeris of 
PSR J0030+0451 from I Abdo eraD ((20091) . 

A s seen from previous X-ray timing studies of this pul- 
sar (jBecker et al.ll2000l : [Becker fc Aschenbachll2002[ ). the 
X-ray pulse proflle of PSR J0030-I-0451 is characterized 
by two peaks that are signiflcantly broader than their ra- 
dio counterparts (Fig. 5). The substantial improvement 
(factor of ~4.5) in photon statistics compared to previ- 
ous observations reveals new information regarding the 
X-ray pulse morphology. For instance, it is now appar- 
ent that the primary pulse is significantly stronger than 
the secondary, with a ~13% larger peak flux in the 0.3-2 
keV range. For thermal polar cap emission this difference 
can be naturally produced by a geometric configuration 
in which one of the polar caps makes a closer approach 
to the line of sight to the observer. There is marginal 
evidence for a larger relative difference between the two 
pulses for 0.7-2 keV compared to 0.3-0.7 keV, although 
statistically they are consistent with being identical. 

Based on the lightcurve model discussed in §6, the 
peak-to-peak separation in phase between the two pulses 
is found to be ^0.55 as measured from the stronger to 
the weaker pulse. It is also evident that the two minima 
between the pulses are uneven, with a fractional differ- 
ence in the minimum fiux of ~22% (for 0.3-2 keV). For a 
centered dipole (i.e. antipodal hot spots), Doppler boost- 
ing and aberration, induced by the rapid stellar rotation, 
can be sufficient to account for the apparent assymetry 
for larger stellar radii (e.g. >12 km for M = 1.4 Mq). 
On the other hand, for more compact stars (e.g. <12 
km for M = 1.4 M©) it would be necessary to invoke 
an off-center dipole in order to reproduce the observed 
asymmetry in the minima as well as the phase separation 
of the two pulses. 

^ Available at |http : / / www. atnf.csiro.au/research/ pulsar / tempo / 
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Fig. 5.— XMM-Newton EPIC pn X-ray pulse profiles of PSR 
J0030-I-0451 for different energy bands. The top panel shows the 
template radio pulse profile at 1420 MHz. The choice of phase zero 
and the alignment of the X-ray and radio profiles are arbitrary. 
Two rotational cycles are shown for clarity. 

To estimate the pulsed fraction of the folded, binned, 
and background-subtracted lightcurves we have used 
two approaches: (i) from the pulse profile fits dis- 
cussed in the following section and (ii) by computing 
th e root-mean-squared pulsed fraction (see, e.g., Eq. 2 
in iDhillon et al.l 120091 for a definition). The former 
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Fig. 6.— XMM-Newton EPIC pn pulse profiles of PSR 
J0030-I-0451 in the 0.3-2 keV {top) and 0.3-0.7 keV and 0.7-2 keV 
{bottom) ranges fitted with a model of a rotating neutron star with 
two-temperature H atmosphere polar caps. The DC level (dashed) 
and uncertainties (dotted) as determined from the fits are shown 
in the top panel. The background level for the 0.3-2 keV band is 
314 counts per bin for 25 phase bins. See text for best fit model 
parameters. 



approach yields 67%+^^° (0.3-2 keV), 71%+^^° (0.3-0.7 
keV), and 63%+^| (0.7-2 keV), whereas using the latter 
method we obtain 64%+^|, 64%+^^°, and 63%+^| for the 
same energy bands, respectively. Although these pulsed 
fractions appear high, they are fully consistent with a 
purely thermal origin of the observed X-rays (as sug- 
gested by the spectroscopic analysis in §4 and as demon- 
strated in §6) if one considers a light-element atmosphere 
at the neutron star surface. Note that the pulse profiles 
in the 0.3-0.7 and 0.7-2 keV energy bands are consistent 
with having the same pulsed fraction. 

Even though the photon statistics are quite limited 
above 2 keV and the profile is consistent with a constant 
flux, the weak modulations in the 2-6 keV band (bot- 
tom panel of Fig. 5) are suggestive of pulsations similar 
to those seen at lower energies. Unfortunately, as a con- 
sequence of the relatively high background level of the 
EPIC pn data, no useful timing information is present 
above '-^3 keV, where the hard spectral tail begins to 
dominate the emission. Thus, we are unable to place 
any limits on the pulsed fraction for 3-10 keV, which 
could potentially offer clues regarding the true nature of 
the hard spectral tail. 

6. CONSTRAINTS ON THE NEUTRON STAR 
COMPACTNESS 



X-rays from PSR J0030+0451 
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The spectrum and pulse profile of PSR J0030+0451 
are indicative of a thermal origin of the observed X- 
rays. This is of particular importance since the surface 
thermal emission can potentially provide useful informa- 
tion regarding the neutron star compactness through a 
measurement of the mass-to-radius rat io (M/R). For 
this p urpose we employ the model from iBogdanov et al.1 
()2007f ). which considers a relativistic rotating compact 
star with two identical, circular X-ray emitting hot spots, 
each (presumably) corresponding to one of the magnetic 
polar caps. The model also incorporates a non-rotating 
Schwarzschild metric as a description of the space-time 
in the vicinity of the star and includes Doppler boost- 
ing and propagation time delays. This relatively simple 
for malism is remarkably accurate for spin periods P ^ 3 
ms ()Cadeau et all [2007t [Morsink et a l. 2007). The NS 
surface is assumed to be covered by a non-magnetic, 
optically-thick H atmosphere. 

Before applying the model to the X-ray pulse profiles, 
it was first convolved with the appropriate instrument 
response, while taking into account the encircled energy 
fraction and the sky and detector background. We fit- 
ted the pulse shape by considering the following param- 
eters: the two temperatures and effective radii of each 
hot spot (Ti, T2, Ri, and R2) the two angles a and 
C, the stellar radius R, the offsets of the secondary hot 
spot from the antipodal position (Aa and A(f)), and the 
phase (j) of the pulse peak. The hydrogen column density 
along the line of sight was fixed at A^h = 2 x 10^° cm~^, 
while the distance was set to D = 300 pc (jLommen et al.l 
[2OO6). Varying A^h over the range (1 - 3) x 10^° cm-^, 
does not result in appreciable changes in the most inter- 
esting parameters, namely M/R, a, and The same 
holds true for the parallax distance uncertainty (~30%), 
which mostly affects the flux normalization, resulting in 
a larger uncertainty in the emission area. As M/R, a, 
and C do not determine the flux normalization, the best 
fit values of these parameters are virtually insensitive to 
the distance uncertainty. Unless noted otherwise, in our 
analysis we assumed a fixed mass of M = 1.4 Mq and 
allow R to vary. The fit was performed simultaneously 
in two energy bands (0.3-0.7 keV, with 0.42-0.5 keV ex- 
cluded as in the spectral fits, and 0.7-2 keV) in order 
to weaken any degeneracies between spectral parameters 
(temperatures and radii) from the geometric and rela- 
tivistic parameters, which have no energy dependence. 
In order to account for any possible presence of a non- 
thermal component, we included the expected number 
of photons per bin for the best fit powerlaw spectrum 
into the uncertainties (assuming an unpulsed or weakly 
pulsed contribution). The best fit values for all free pa- 
rameters were found by searching the hyperspace. To 
ensure that the obtained values correspond to the abso- 
lute minima of the parameter space, the fit was repeated 
for 100 different combinations of initial parameter values. 

The resulting best fit pulse profiles for a H atmosphere 
model are presented in Figure 6. Unlike blackbody emis- 
sion, the atmosphere model is able to reproduce the pulse 
shape and large amplitude remarkably well. The stel- 
lar radius is constrained to be i? > 10.4 km or, more 
generally, R/Rs > 2.5 (99.9% confidence) for all com- 
binations of the other free parameters. In the case of 
a blackbody model, good fits to the X-ray pulse pro- 
file require implausably large stellar radii (j^20 km) for 
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Fig. 7. — The mass-radius plane for neutron stars showing various 
theoretical model tracks. The regions ruled out by the X-ray data 
for PSR J0030+0451 at 95% and 99.9% confidence are shown in 
blue and green, respectively. 

a range of assumed massess (M = 1 — 2 Mq). This is 
in ag reement with the resu lts obtained for PSR J0437- 
4715 (IBogdanov et al.|[200l . PSR J2124 -3358, and using 
the a rchival data for PSR J0030-f 0451 (|Bogdanov et al.1 
[2OOI . 

Figure 7 shows the M — R plane with a sample 
of predicted model track for different NS EOS from 
iLattimer fc PrakashI ()2001[ ) and the constraints obtained 
for PSR J0030+0451, with the shaded regions corre- 
sponding to regions excluded by the limits on M/R. Al- 
though the results are consistent with all EOS if M is 
unconstrained, it is still possible to draw some impor- 
tant conclusions. In particular, there is observational 
evidence that several recycled p ulsars are niore ma s- 
sive than the canonical 1.4 MO (iFreire et al.l l2008allbl : 
IVerbeist et aIll20Q8l: IChampion et al.ll2008D . as expected 
from standard MSP formation theory, which involves ac- 
cretion of an appreciable amount of material. If this is 
true for isolated MSPs such as PSR J0030-f 0451 as well, 
then for M > 1.45 the M/R constraint is inconsistent 
with some kaon condensate EOS (e.g. GSl for which the 
maximum mass consistent with our limits is ^1.3 Mq), 
and most quark star EOS (e.g., SQM1-SQM3 in Fig. 7). 
For M > 1.6 Mq, the limit of R/Rs > 2.5 rules out aU 
but the stiffest equations of state (such as MS0-MS2). 
Conversely, if this pulsar is in fact a quark star, then 
the constraints obtained herein would require it to be 
M < 1.45 Mq. 

The fits to the pulse profile also provide limits on the 
global magnetic field and viewing geometry of the pul- 
sar. Specifically, the pulsar obliquity is constrainted to 
be in the range a — 38° — 90°, implying a substantial 
misalignment between the spin and magnetic axes, while 
the angle between the line of sight and the spin axis is 
C = 47° — 90°. The relatively broad range of allowed 
values arises due to the strong covariance of these an- 
gles with the offset in latitude of the secondary polar 
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cap Aa, which is weakly constrained to be in the range 
—42° — 52°. As defined, a positive value of Aa indicates 
a southward offset from the antipodal location. The dis- 
placement in the longitudinal direction is constrainted to 
be A0 = — 8° — 46°, consistent with both a centered and 
offset dipole. For this parameter, a positive value corre- 
sponds to an offset in the direction of rotation, meaning 
that the hot spot leads the antipodal position. 

7. CONCLUSION 

In this paper, we have presented deep XMM- 
Newton observations of the nearby recycled pulsar PSR 
J0030-I-0451. For the first time we are able to fully char- 
acterize the emission spectrum of this isolated pulsar in 
the 0.1-10 keV band. It exhibits striking qualitative sim- 
ilarities to that of PSR J0437-4715, suggesting that other 
typic al MSPs such as the nearby J2124-3358 and J1024- 
0719 lZavlin"2006'), as well as the majority of MSPs in 47 
Tuc ([Bogdanov et al. 2006a), whether binary or isolated, 
may exhibit the same characteristic spectral shape over 
this energy range. 

Although the spectral continuum of J0030-I-0451 can 
be adequately described by a variety of plausible mod- 
els, in all instances the bulk of detected photons in the 
soft X-ray range (0.1-2 keV) appear to be of thermal 
origin. Our H atmosphere polar cap model is in fair 
agreement with the observed continuum and pulsed emis- 
sion from PSR J0030+ 0451. As with PSR J0437-4715 
(|Bogdanov et al.l [20071 ) . the relatively large pulsed frac- 
tions require the existence of a light-element atmosphere 
on the stellar surface and cannot be reproduced by a 
blackbody model for realistic NS radii. By modeling the 
thermal pulse shape we are able to place a limit on the 
allowed stellar radii of i? > 10.4 km (99.9% confidence) 
assuming M = 1.4 Mq, a substantial improvement over 
the constraints obtained from the shorter arch i val ob - 
servations. As demonstrated in iBogdanov et all (|2008f ) . 
even deeper X-ray observations of this and other MSPs 
should lead to stricter limits on M/R, which in turn may 



firmly rule out entire families of neutron star EOS. This 
is especially true for the growing number of binary MSPs 
for which complementary mass measurements from radio 
pulse timing are becoming available. 

With future observations, it is important to uncover 
the true nature of the hard emission (> 3 keV), through 
a combination of deep energy-resolved timing and phase- 
resolved spectroscopy at energies above 3 keV. Better 
characterization of the optical and 7-ray spectra of PSR 
J0030-I-0451 could potentially provide additional infor- 
mation regarding the origin of this emission. These 
endeavors have important implications for establishing 
MSP X-ray timing as a vible means for tight constraints 
on the neutron star EOS. The potential utility of MSP 
as precision probes of neutron star structure makes them 
particularly important targets for upcoming X-ray mis- 
sion such as the International X-ray Observatory. The 
great increase in sensitivity of such a facility would per- 
mit detailed observations of a larger sample of MSPs than 
currently possible, which could ultimately lead to defini- 
tive insight into the poorly understood properties of the 
neutron star interior. 
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